A remark on Krein's resolvent formula and boundary 
conditions 
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Krein's resolvent formula 1 1 1 is a powerful tool in the spectral analysis of self-adjoint 
extensions, which found numerous applications in many areas of mathematics and physics, 
including the study of exactly solvable models in quantum physics 0|3]|4). For the use of 
this formula in the traditional way one needs a kind of preliminary construction, like finding a 
maximal common part of two extensions, see [2 Appendix A]. While this is enough for many 
applications, including models with point interactions, there is a number of problems like the 
study of quantum graphs or more general hybrid structures, where self-adjoint extensions are 
suitable described by more complicated boundary conditions, see |5]|6]0, and it is necessary 
to modify Krein's resolvent formula to take into account these new needs. This can be done 
if one either modifies the coordinates in which the boundary data are calculated (6| |8] or 
considers boundary conditions given in a non-operator way using linear relations (9] 1101 . 
On the other hand, a more attractive idea is to have a resolvent formula taking directly 
the boundary conditions into account, without changing the coordinates. We describe the 
realization of this idea in the present note. 

Let S be a closed densely defined symmetric operator with the deficiency indices (n,n), 
< n < °°, acting in a certain Hilbert space Jf . One says that a triple (V,ri,r2), where 
V = C" and F\ and Y2 are linear maps from the domain domS* of the adjoint of S to V, 
is a boundary value space for S if (0,5*1//) — (S*<j),y) = (ri0,r2V) — (r20,r2V) for any 
0, y £ domS* and the map (Ti ,1^) : domS* — > V V is surjective. A boundary value space 
always exists [9 Theorem 3.1.5]. All self-adjoint extensions of S are restrictions of S* to 
functions £ domS* satisfying AT\§ = where the matrices A and B must obey the 

following two properties: 



AB* = BA* AB* is self-adjoint, 

the n x 2« matrix (A\B) has maximal rank n. 



(1) 
(2) 
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We denote such an extension of S by H A ' B . Our aim is to write a formula for the resolvent 
R A,B (z) — (H A - B — z)^ 1 in terms of these two matrices A and B. 

We will need some notions from the theory of linear relations. Let V = C. Any linear 
subspace A of V © V will be called a linear relation on V. By the domain of A we mean the 
set domA = {x G V : By G V with (jt,y) e A)}. A linear relation A -1 = {(x,y) : (y,x) € A} 
is called inverse to A. For a G C we put aA = {(x, ay) : (jc,y) € A}. For two linear relations 
A', A" C V © V one can define the sum A' + A" = {(*,/ +/'), {x,y') G A', (*,/') e A"}; 
clearly, one has dom(A' + A") = domA' n domA". The graph of any linear operator L acting 
in V is a linear relation, which we denote by grL. Clearly, if L is an invertible operator, then 
grL -1 = (grL) -1 . For arbitrary linear operators L',L" one has grL' +grL" = gr(L' +L"). 
Therefore, the set of linear operators is naturally imbedded into the set of linear relations. 

Denote by J an operator acting in V ©V by the rule J(xi,x 2 ) — (x 2 ,—Xi), X\,x 2 \inV. 
For a linear relation A on V the relation A* = /A x is called adjoint to A; A is called 
symmetric if A C A* and is called self-adjoint if A = A*. The graph of a linear operator 
L in V is symmetric (respectively, self-adjoint), iff its graph is a symmetric (respectively, self- 
adjoint) linear relation. In other words, a self-adjoint linear relation (abbreviated as s.a.l.r) 
is a symmetric linear relation of dimension n. Let A, B be n x n matrices. We introduce the 
notation 

A AM = {(x u x 2 ) GVffiV, Ax l =Bx 2 }. 

A criterion for A A B to be self-adjoint was proven in Q: A linear relation A A ' B is self- 
adjoint iff A and B satisfy {1} and (0. It is important to emphasize that any s.a.l.r. A 
can be defined by this construction, more precisely, there exists a unitary operator U such 
that A = A'( 1+c/ ) ,1 ~ £/ (5|. This shows that there is a bijection between s.a.l.r.'s and unitary 
operators; nevertheless, we consider parametrization by the two matrices of coefficients as a 
more natural way to present s.a.l.r. 

The language of linear relations is widely used in the theory of self-adjoint extensions 
of symmetric operators ll3ll8l ll II . Let us return to the symmetric operator S and its boundary 
value space (V,ri,r2). It is a well known fact that there is a bijection between all self- 
adjoint extensions of S and s.a.l.r's on V. A self-adjoint extension H A corresponding to a 
s.a.l.r. A is a restriction ofS* to elements <j) G domS* satisfying abstract boundary conditions 
(ri0,r20) G A 1 9 Theorem 3.1.6]. To carry out the spectral analysis of the operators H A it 
is useful to know their resolvents, which are provided by the famous Krein's formula | lI llOI . 
To write this formula we need some additional constructions. For z G C \ R, let J/ z denote the 
corresponding deficiency subspace for S, i.e. JV Z — ker(5* — z). The restriction of T\ and Y 2 
onto Jf z are invertible linear maps from ,yK to V. Put y z — (T\ \jr\ 1 and Q(z) = T 2 y z ; these 
maps form holomorphic families from C\M to the spaces _£? (V,Jrf?) and Jz? (V, V) of bounded 
operators from V — > Jrf? and from V to V respectively. Denote by H° a self-adjoint extension 
of S given by the boundary conditions =0. i.e. for A = 1 and B = 0, then the maps y z 
and Q(z) have analytic continuations to the resolvent set res// , and for all z, £ G res// one 
has CHI 

Q(z)-Q*(Q = (z-Qr^Yz- (3) 

The maps y, and Q(z) are called the T-field and the i?-function for the pair (S,H°), 
respectively. 

For all z G res// nres// A we consider the linear relation gr<2(z) — A. While this set is, 
generally speaking, not the graph of an operator, the inverse linear relation (gr Q(z) — A) 1 
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is the graph of a certain linear operator C A (z), so that the resolvent R A (z) = (H A — z) 1 is 
expressed through the resolvent R°(z) = (H°—z)~ l by Krein's formula H 1 01 Proposition 21 

R A (z)=R°(z)-YzC A (z)Y,*. (4) 

The calculation of C (z) is a rather difficult technical problem, as it involves "generalized" 
operations with linear relations. Such difficulties do not arise if A is a graph of a certain linear 
operator L; the corresponding boundary conditions can be presented by 

r 2 0=Lr^, (5) 

and such extensions are called disjoint with respect to H° because they satisfy the equality 
dom// A n domi/ = domS Q] (the operator S is then called the maximal part of H° and 
H A ). Then the subspace giQ(z) — A is the graph of the invertible operator Q(z) — L, and 
C A (z) = (Q(z)-L)~. Actually, for a given self-adjoint extension H one can find a boundary 
value space (which is, of course, not unique) such that H corresponds to the boundary 
conditions (0 with a suitable L (6] [8). But finding such boundary value space involves a 
lot of other problems, in particular, the operators Q(z) and y, must be changed 1 10 8 12 1. On 
the other hand, any boundary conditions can be represented with the help of two matrices by 

AT 1 <j>=Br 2 <l> o (rK/>,r 2 </>)eA A ' B (6) 

with A and B satisfying Q and (|2j- Our aim is to show that the resolvent formula admits a 
simple form in terms of these two boundary matrices. Here is the main result of our note. 

Proposition (Modified Krein's resolvent formula). Let H AB be the self-adjoint extension 
of S corresponding to the boundary conditions with A, B satisfying Q and (|2}, and 
z £ resi/ Dresi/' 4 ' B , then 

(a) the matrices Q(z)B* — A* and BQ{z) —A are non-degenerate, 

(b) the resolvent R A B (z) = (H A B - z) _1 is connected with R°(z) by 

R A ' B (z)=R (z)-y z B*(Q(z)B*-A*y 1 yt, (7) 
or 

R A > B (z)=R°(z)-y z (BQ(z)-A)- l B Y i (8) 

Remark. The formulas and © are equivalent. To obtain them from each other one should 
replace z by z and take in both sides adjoint operators taking into account the resolvent 
property R(z) — R*(z) and the equality Q(z) — Q*(z) which follows from 0. 
First we prove some simple properties of the matrices A and B: 

Lemma. Let A, B satisfy Q and (0, then 

(a) kerA*nkerfi* =0, 

(b) A AM = {(B*x,A*x), x e V}. 

Proof of Lemma, (a) We first remark that the condition Q is equivalent to ranA + ranZ? = V. 
Then kerA* nkerB* = (ranA) x n (ranB) x = (ranA + ranB) x = V x = 0. 

(b) Eq. says that {(B*x,A*x), ieV}c A AB . At the same time, it follows from (a) 
that the linear subspace on the left-hand side has dimension n, which coincides with the 
dimension of A A - B . Therefore, these two linear subspaces coincide. □ 
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Proof of Proposition, (a) The matrices in question are adjoint to each other, therefore, it 
suffices to prove that one of them is non-degenerate. 

Consider first the case zeC\R. As follows from (|3jl, for such z the matrix Imz • Im Q(z) 
is positive definite (Here ImQ = {Q- Q*)/(2i).) Assume that det (Q(z)B* -A*) = 0, then 
there exists a non-zero x E V with 

{Q{z)B*-A*)x = 0. (9) 

If B*x = 0, then, due to item (a) of Lemma, we would have A*x ^ 0, and l|9} would be 
impossible. Therefore, B*x ^ 0. Taking the scalar product of B*x with both sides of l|9} we 
get (Q(z)B*x,B*x) = (x,AB*x). The number on the left-hand side has non-zero imaginary 
part, while the number on the right-hand side is real due to (0. This contradiction proves the 
requested non-degeneracy for non-real z- 

NowletzeMnresi/ anddet(2(z)B*-A*) =0. Let us show that zE spec H AB . Eq. (0 
says that Q(z) is now self-adjoint. As the matrix Q(z)B* — A* in non-invertible, there is a non- 
zero x E ran(g(z)B* — A*) X = ker(fig(z) — A). By definition, the element <j> = y z x is an 
eigenvector of S* corresponding to the eigenvalue z- Let us show that <j) E domH AB (and then 
z is an eigenvalue ofH AB ). In fact, one has r\<f> = T\y z x = TiT^ x =x, r 2 = = Q(z)x, 
and, therefore, AT\x — BTix = — (BQ(z) — A)x = 0, which means that (j) £ domH AB . 

(b) Taking into account the remark after Proposition it is enough to prove only 0. 
Actually, we only must show that for z E resH n resH AB the linear relation inverse to 
grQ(z) — A A B is the graph of B* (Q(z)B* — A*) \ Taking into account item (b) of Lemma, 
we conclude that domg(z) n domA A ' B = ranfi* and grg(z) - A A - B = ^{B*u,(Q(z)B* - 

A*)u),uev} = {(B*(e(z)B*-A*)-Vw),we v}. □ 

Let us consider some example from the point of view of the resolvent formula. 

Example (Graph with a single vertex). Let = ®" =l Jfj with Jtfj = _Sf 2 (R^), where 
each rV' is a copy of the positive half-line [0,+°°) . We will write the elements of 
€ Jif in the vector form, <p(x) = ($j(xj)), <pj E J&j, Xj E . By S we denote an 
operator which acts on each Jffj as —d 2 /dx 2 with the domain domS = {0 = (0 ; ), (j>j £ 

W 2 ' 2 (R^ -) ), <j)j(0) = <j)'j(0) =0,j=l,.. .,«}. Then S* is the direct sum ® n j=l d 2 /dxj with the 
domain Q)" =l W 2,2 (M.^ . An integration by parts show that as a boundary value space one can 

take (C",r l5 r 2 ) within =i\ = (-4>j(o)) = -f (o), r 2 =r 2 (to) = (^(o)) = 0(0); 

details may be found in j5j|6][7). Then any self-adjoint extension of S involves the boundary 
conditions A<j>'(0) +B<p(0) = with suitable A and B, which describes the coupling of n 
half-lines at the origin. The operator H° defined by 0'(O) = is the direct sum of Neumann 
Laplacians, and its Green function (the resolvent integral kernel) G°(x,y;z) is given by 

G°(x,y;z) = — ^= diag (e~^~^^i~ y ^ + e -^ x J +y ^), 
2y — z 

where the contiuous branch of the square root is chosen by the rule Re y/z > for z ^ (— °°, 0] , 

and x = {xj), y — (yj), xj,yj E M.+ . Then the elements g{ = G°(-,0;z)ej = e~*S-~ zx / 's/^ze j, 
(ej) is the standard basis in C", j = 1, . . . , n, form a basis in jY z = ker(5* — z), and the map 
y z defined by C" 9 v = (vy) i— > (vjgl) is the corresponding T-field, because T\ y z v — v for any 
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v S C", and the ^S-function has the simple form Q(z) — 1 / y/^zE„, where £„ is the unit matrix 
of order n. Therfore, the resolvent for H A B takes the form 

R A ' B (z)=R°(z)- t Cjt{z){g k z M 

with C(z) = (BQ(z) -A)' 1 B = ^/^z(B - ^/^zA)- l B. A similar formula was obtained in 0. 

Example (Point interactions with mixed boundary conditions). Here we consider the well- 
known example of point interactions in three dimension |2 1. Let Y = (y\, . . . ,y„) C R 3 . Denote 
by S the Laplacian with the domain domS = {0 g W 2 ' 2 (R 3 ), 0(F) = 0}. This operator 
has deficiency indices (n,n), and the adjoint operator is the Laplacian with the domain 
domS* = W 2;2 (R 3 \ Y). Each function <j) e domS* has the following asymptotics: 

${x) = ^-j-^}+<pj + o(l), x^yj, &C 1 , j = l,...,n, 

and the vectors — and = (0?) can be considered as boundary values of 0, 
see I3ll6l ll3l for details. The operator// is just the free Laplacian in R 3 . Denote by G°(x,y;z) 
its Green function, G°(x,y;z) — e~ v/Z ^ x ~- v '/(4ft|x — y\), then the functions g{ = G°(-,y/;z), 
j = 1 , . . . , n, form a basis in ker(5'* — z), and the map y z : C"9v = (v ^ ) i— > £ v /g^ is the T-field 
(it is easy to check that Y\ y z y — v), and the J?-function is given by the nx n matrix 

QAZ) = G^y k ; Z ) = J / *, QjjV = 

Therefore, the resolvent of the operator H A B given by the boundary conditions A0 1 = B(j> 2 
can be defined by its intergal kernel 

G A ^(x,y-z)=G°(x,y-z)~ £ C ;1 (z)G (x,y ; -;z)G°( W ,y;z) 

with C(z) = — -A) B. We remark that the class of interactions described by this 

formula is wider that the one studied in 0. Some properties of H AB in its dependence on A 
and B were studied recently in 1131 . 
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